arXiv: 1506.01196v2 [astro-ph.HE] 18 Mar 2016 


Cosmic-Ray Hardenings in the Light of AMS-02 

Yutaka Ohira/’B Norita Kawanaka,^ and Kunihito loka^’^ 

^Department of Physics and Mathematics, Aoyama Gakuin University, 5-10-1 Fuchinobe, Sagamihara 252-5258, Japan 

^Department of Astronomy, Graduate School of Science, 

University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan 
^ Theory Genter, Institute of Particle and Nuclear Studies, KEK, Tsukuba 305-0801, Japan 
* Department of Particle and Nuclear Physics, SOKENDAI (the 
Graduate University for Advanced Studies), Tsukuba 305-0801, Japan 

Recent precise observations of cosmic rays (CRs) by AMS-02 experiment clearly show (1) harder 
spectra of helium and carbon compared to protons by oc and (2) concave breaks in proton 

and helinm spectra at a rigidity R ~ 300 GV. In particular the helium and carbon spectra are 
exactly similar, pointing to the same acceleration site. We examine possible interpretations of these 
features and identify a chemically enriched region, that is, superbubbles as the most probable origin 
of Galactic CRs in high rigidity i? > 30 GV. The similar spectra of CR carbon and helium further 
suggest that the CRs with i? > 30 GV originate from the supernova ejecta in the superbubble core, 
mixed with comparable or less amount of interstellar medium. We predict similar spectra for heavy 
nuclei. 

PACS numbers: 98.70.Sa, 98.38.Mz, 96.50.sb 


I. INTRODUCTION 

Supernova remnants (SNRs) are thought to be the ori¬ 
gin of Galactic cosmic rays (CRs). Yet, we have not un¬ 
derstood what type of supernova dominantly produces 
Galactic CRs in each energy, isolated supernovae 
or supernovae in superbubbles [l|. There are many sug¬ 
gestions and interesting discussions in previous studies 
about CR composition^ that are based on CR data in 
the rigidity R < 1 GV 0 . Now, we consider high energy 
CRs with R 1 GV. If the spectral index depends on 
CR species, the CR composition should depend on the 
rigidity, R. 

Recently accurate measurements by AMS-02 experi¬ 
ment grea^ improve previous results by PAMELA 
[1], CREAM dl, AMS and ATIC-2. Two kinds of ”CR 
hardening” are observed by AMS-02: 

1. Spectra of CR helium and carbon are harder than 
that of protons, and the spectral index of CR car¬ 
bon is almost the same as that of CR helium (see 
Figure [T]). 

2. The spectra of CR protons and helium have an up¬ 
turn break at a rigidity R ^ 300 GV (see Figure[2]). 
The CR carbon possibly has the same break al¬ 
though the data is not sufficient. 

The precision measurements of carbon up to i? 1 TV 
are newly provided by AMS-02. From Figure [TJ we find 
that the ratios, He/p and C/p, are well approximated 
by the power law as oc Although the spallation 

and solar modulation reduce helium and carbon at low 
energy, they distort the power law by less than ^ 20% 
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FIG. 1. The helium to proton ratio (He/p; filled circle) and 
the carbon to proton ratio (C/p; filled sguare) as a function of 
rigidity R, that is provided by AMS-02 [^. These ratios are 
not constant but clearly rise as In addition, these ra¬ 

tios are appreciably larger than the solar abundance {dashed 
lines). On the other hand, the carbon to helium ratio is con¬ 
stant. 

in the energy range of 10-10^ GV where the effect 
of spallation may be estimated from the boron-to-carbon 
ratio (B/C). Remarkably, the CR carbon has the same 
spectrum as the CR helium, which strongly suggests that 
both carbon and helium are accelerated at the same place 
(see Section IIVDI for detail). Furthermore, the abun¬ 
dances of the CR carbon and helium are larger than the 
solar abundance Q- Especially for the CR carbon, the 
abundance is over 10 times larger. 

Before AMS-02, several models were proposed to ex¬ 
plain the helium hardening. The new AMS-02 data, es- 
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FIG. 2. Spectra of CR protons (red), helium (green), and 
carbon (blue) fluxes as a function of rigidity. The data are 
given by AMS-02 and CREAM [^. These are fitted by 
theoretical curves based on the superbubble model in Section 
2 (sold lines). 


pecially of carbon, can narrow down these models and 
the CR sources. In this paper, we discuss which model is 
better to explain the new data and conclude that the su¬ 
perbubble model is most likely (see Table |T|. Therefore, 
we first briefly show that our superbubble model Q can 
explain the new data of AMS-02 (Section [HI. Then, we 
discuss which model is better to explain several features 
of observed CR spectra, i.e., the overabundance of CR 
carbon (Section IIIII) , the different spectra between CR 
protons and CR helium/carbon (Section |TVl Table HI, 
and the spectral breaks at i? ~ 300 GV (Section IVT) . 
Section |VT] is devoted to summary and discussion. 

Note that AMS-02 also provides antimatter data. Here 
we assume that the antiprotons and positrons are pro¬ 
duced by other sources, such as nearby SNRs If 
nearby SNRs interact with dense molecular clouds, sec¬ 
ondary antiprotons and positrons originated from the 
SNRs can dominate the Galactic component, while GR 
protons from them fall below the Galactic component. 


II. SUPERBUBBLE MODEL 

According to the standard theory of Galactic GRs, GRs 
are accelerated by SNRs. The CRs escape from SNRs to 
the interstellar medium (ISM), propagate in our Galaxy, 
and finally escape from our Galaxy. The steady-state 
spectra of GR protons, helium and carbon are given by 

OC R ’ , Sobs,2 — T Ti : (1) 

where the subscript i = p, He, G represents particle 
species and 7 ^ is the rigidity dependence of the diffusion 


coefficient in our Galaxy, D oc W'. Sinj,i is the spectral 
index of GRs injected into the ISM. 

High energy GRs are thought to be accelerated in SNRs 
and escape from them earlier than low energy CRs be¬ 
cause the shock velocity is fast and the magnetic field 
is strongly amplified in the early phase. On the other 
hand, low energy CRs are mainly accelerated in SNRs 
at later times because the mass swept by the shock of 
SNRs increases with time and CRs injected early lose 
their energy by the adiabatic expansion before they es¬ 
cape. Therefore, the spectra of CRs inside SNRs are 
given by 


dNi o 

li “ 


( 2 ) 


and the time evolution of the maximum rigidity is given 
by 


Rmax,2 OC t . (3) 

ai and fii are parameters to describe the evolution of the 
maximum rigidity and the number of CRs, respectively, 
and Si is the spectral index of CRs inside SNRs. By 
taking account of the rigidity dependence of CR escape 
from sources [l^, the index, Sinj,i, is given by 





( 4 ) 


If all CRs originate from the same source, the physics 
parametrized by Si, ai and 7 ^ depends only on particle 
rigidity. Then, the inferred rigidity spectrum should be 
universal. However, the number evolution parameter, /3i, 
could depend on the particle species if the composition 
distribution is not uniform because Pi depends on the 
density distribution around SNRs. If so, the steady-state 
spectra of CRs (Equation [T]) depend on the CR species. 

In addition, taking account of the rigidity dependence 
of CR escape from sources 0 , the spectral index, Si, in 
Equation (j4]) may depend on the rigidity because it de¬ 
pends on the Mach number of SNR shocks, M. Accord¬ 
ing to the standard diffusive shock acceleration (DSA) 
theory m , the index, Si, is given by 
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- 1 


( 5 ) 


High energy CRs escape from the CR sources when M 
is large, and Si = 2, while low energy CRs escape later 
when M becomes low, and Si > 2. Then we expect the 
spectral breaks for all CR species at the same rigidity. 
It should be noted that if one of Si,ai,Pi and 7 ^ has a 
rigidity dependence, the steady-state CR spectra have a 
spectral break (see Section IVT) . 

To explain both the helium harding and the concave 
breaks at i? ^ 300 GV, Ref. Q considered a hot gas with 
temperature of T = 10® K and the density profile of 
helium, 

riHe oc r 


( 6 ) 
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TABLE I. Score Sheet for Models of Different Cosmic-ray Spectra 


Model Mechanism Score Comment Reference 


Propagation 

Spallation 

C 

Inconsistent with data of B/C 

US 

Different sources 

Wind 

C 

Requires suppression of He (p) acceleration in the ISM (winds) 

LIS 


Reverse shock 

B 

Requires suppression of He acceleration in the ISM 

18J 

Injection 

Injection 

C 

Inconsistent with C/He=const. 

19J 

Inhomogeneous environment 

Ionization 

C 

Inconsistent with C/He=const. 

20] 


Superbubble 

A 

Consistent with observations, and CRs originate from ejecta. 

[8J 


being different from that of protons, where r is the ra¬ 
dius of a system. Such plasma is realized in superbubbles 
(see Sections llVDI) . Figure [2] shows that our superbub¬ 
ble model is in excellent agreement with the new data 
of AMS-02 by using parameters a = 6.5,7 = 0.44, (5 = 
0.52, EHe = 0.4, ec = 0.013, and T = 10® K (see Ref.[| 
for detail). The similar spectra of helium and carbon 
suggest similar density profiles of these elements. 


III. CARBON OVERABUNDANCE 

There is a standard picture for the carbon overabun¬ 
dance based on the past observations of CR compositions 
[3 [11. The observational fact is that the abundance 
of refractory elements relative to solar is systematically 
larger than that of volatile ones. The interpretation of 
this behavior is that the refractory elements condense 
into dust grains and the dust grains are injected into the 
acceleration more efficiently than interstellar gas-phase 
ions because of their enormous gyroradii. The sput¬ 
tering of accelerated dust yield seed ions and they are 
accelerated to the observed CR nuclei with overabun¬ 
dance. On the other hand, among the volatile elements, 
the abundance enhancements relative to solar increase 
with mass, whereas the refractory component shows lit¬ 
tle mass dependence. The mass dependence of volatile 
elements may reflect the scattering cross section between 
fast grains and the ambient gas [l4l| or a mass-to-charge 
dependence of the injection efficiency via the gyroradius 
filtering. Thus the volatile elements are presumably in¬ 
jected directly from the interstellar gas. 

CR carbon falls between the refractory and volatile 
elements. Therefore it is reasonable that the CR carbon 
is injected from dust grains as well as gas-phase ions. 
This provides an important clue to the origin of CRs as 
discussed in the next section. 


IV. ORIGIN OF DIFFERENT SPECTRA 

It is a mystery that the spectral index of protons 
is different from that of helium and carbon because 
the same rigidities give the same gyroradii and hence 
the usual acceleration and propagation are conmosition 
blind. Since the helium hardening was observed [1,S,[3|, 
several models have been proposed for explaining the 


spectral difference between helium and protons. These 
can be broadly divided into four types: propagation 
ISection riV Ajl . different sources (Section |TVB|), injection 
physics 1 Section II V C|) and inhomogeneous environment 
fSection IIVD|) . where the last three are related to the 
nature of the CR accelerators. 

The carbon spectrum precisely observed by AMS-02 
provides a hint to the riddle of the different spectra. The 
spectral index of carbon is almost equal to that of he¬ 
lium and strongly suggests that the acceleration site is 
the same for helium and carbon. In addition the amount 
of injection should keep a constant ratio between helium 
and carbon even if carbon dust coexists with helium gas 
(see Section uni). These requirements severely constraint 
the physical condition of the acceleration site. In the fol¬ 
lowing we inspect the models for the spectral differences 
in the light of the AMS-02 results. The score sheet for 
each model is summarized in Table HI 


A. Propagation 

The CR compositions do change duri ng p ropagation 
through the interactions with the ISM [^. However 
spallation of helium and carbon is not effective under the 
grammage inferred from B/C and the antiproton flux [^. 
Thus, hereafter, we discard this possibility and concen¬ 
trate on the models related with the CR accelerators. 


B. Different sources 

The different spectra may simply reflect different 
sources, i.e., the proton sources may be naively differ¬ 
ent from the helium and carbon sources. For this type 
of models, in general, the following conditions should be 
satisfied: 

1 . A fine-tuning is necessary in the sense that com¬ 
pletely different sources realize the helium to pro¬ 
ton ratio comparable to the solar abundance within 
a factor of two. 

2. In the helium (proton) sources, the proton (helium) 
acceleration should be suppressed below the solar 
abundance. Otherwise the helium to proton ratio 
would deviate from a power law. 
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3. The position of the spectral breaks should coincide 
between protons and helium even if the sources are 
different. This is not trivial since the physical con¬ 
ditions at the acceleration sites are different. Thus 
we may have to call for propagation to produce the 
spectral breaks (see Section 

In the real situation, these conditions seem difficult to 
realize. One representative model is the wind model, in 
which a supernova explodes into the stellar wind and the 
CR helium and carbon are accelerated in the wind, while 
the protons are accelerated from ISM [I3- Helium and 
carbon are abundant in the wind material if the cen¬ 
tral star is the progenitor of type Ib and Ic supernovae, 
respectively. However type H supernovae occur most fre¬ 
quently and the composition of their winds is quite sim¬ 
ilar to that of ISM. Thus, the proton acceleration in the 
wind of type H supernovae should be suppressed below 
the solar abundance by an unknown mechanism. Fur¬ 
thermore, in this model the helium acceleration in the 
ISM should be suppressed below solar abundance by an¬ 
other unknown mechanism. Therefore, it seems hard to 
explain both suppressions simultaneously. 

The other model is the reverse shock model, in which 
the reverse shock running back into the supernova ejecta 
accelerates the helium and carbon [l^. Note that type la 
supernovae do not produce helium and carbon so much. 
However the reverse shock is much less energetic than 
the forward shock except for a very brief period of the 
shock-crossing and accelerated particles may suffer from 
adiabatic losses. Moreover, as with the wind model, the 
helium acceleration should be somehow suppressed below 
solar abundance in the forward shock by an unknown 
mechanism. 


C. Injection physics 

In the injection model, the differences between ele¬ 
ments arise from the injection process in the collisionless 
shock acceleration M- According to the injection model 
of Ref. El, the injection into DSA is more efficient at 
smaller Alfven Mach numbers, M^, and protons have the 
strong dependence on Ma compared with helium. Be¬ 
sides, the spectra of accelerated particles are softer when 
the shock is weaker. Then the time-integrated spectrum 
of CR protons becomes softer than that of helium, as 
observed. 

However, if this is the case, the carbon to helium ratio 
deviates from a constant because the rigidity of carbon, 
especially carbon dust (See Section Ell and partially- 
ionized carbon ions, is different from that of helium. This 
is inconsistent with the AMS-02 observations. Thus the 
injection physics is not likely the cause of the different 
spectra. 

Note that the rigidity of carbon ions are equal to that 
of helium if carbons are fully ionized in high temperature, 
and the carbon overabundance may be also explained if 
CRs originate from carbon rich plasma like superbubbles. 


In this case, however, we should change the standard 
picture for the dichotomy between refractory and volatile 
elements in Section EH 


D. Inhomogeneous environment 

The energy spectra may differ between elements if the 
composition of the ambient medium is not uniform (8l.f20l|. 
Because a stronger shock accelerates higher-energy CRs, 
the outward-decreasing abundance naturally leads to the 
hard spectrum of CR helium. One can easily under¬ 
stand this mechanism in the escape model in which the 
highest-energy CRs escape from the source and shape 
the observed energy spectrum. Higher-energy helium is 
produced in the inner region where the helium is more 
abundant and the shock is stronger. 

The composition gradient may be produced by a gra¬ 
dient of the ionization degree of helium if neutral helium 
atoms are not injected into acceleration [2^. However, 
neutral atoms could be injected into acceleration [El- 
Moreover, the carbon (dust and partially-ionized ions) is 
expected to have a different ionization profile from that 
of helium. Thereby the carbon to helium ratio is not kept 
constant. This is inconsistent with the AMS-02 observa¬ 
tions. 

The chemical inhomogeneity is expected in superbub¬ 
bles. A superbubble is created by multiple supernovae 
and stellar winds in OB associations. The supernova 
ejecta dominates the bubble mass especially in the core 
region [2^. Since the supernova progenitors have too 
small proper motion to escape the bubble, the super¬ 
nova shock may sweep the ejecta-enriched matter. We 
emphasize that the ejecta-enriched region is essential for 
the helium inhomogeneity because the stellar nucleosyn¬ 
thesis can not double the cosmic mean helium abundance 
because of the Schonberg-Chandrasekhar limit This 
is one of the reason why we need the big bang nucleosyn¬ 
thesis. Therefore the superbubble model with concentra¬ 
tion of composition is an inevitable consequence of our 
arguments. 

To reproduce the constant carbon to helium ratio as 
observed by AMS-02 in Figure E the helium and carbon 
should be well mixed in the supernova ejecta such as via 
the Rayleigh-Taylor instability or in the superbubble via 
supernova-driven turbulence. Then we also expect simi¬ 
lar spectra for heavy metals (possibly except for CR iron 
which might be produced by type la supernovae). On the 
other hand, the proton fraction may be different since the 
hydrogen is ejected at completely different spacetime and 
velocities. 

More importantly, we can argue that the old ISM with 
solar abundance should not mix with the fresh super¬ 
nova ejecta so much in the acceleration region. Be¬ 
cause primordial helium and carbon mass fractions are 
Yp = 0.2477 ±0.0029, Zc^p = 0 [2^ while those fractions 
in the present ISM are Iism = 0.275, ^c.ism = 0.0034245 
[ 2 ^, the mean mass ratio of carbon to helium in stellar 
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ejecta is given by 


■^c,ej 




^c,ISM 
biSM — 


0.1254 , 


(7) 


which is much larger than that of ISM, 2'c,isM/hjsM = 
0.01245. In order to produce a hard spectrum of CR he¬ 
lium compared to that of CR protons, our model needs a 
helium rich region at the center of superbubbles. There¬ 
fore, the ejecta should dominate over ISM at the center 
of superbubbles and the carbon to helium ratio, Z^jY, 
should be close to Zc^ej/^ej- If the mass fraction of ISM 
is larger than that of stellar ejecta in the outer region of 
superbubbles, Z^jY should be close to ^c,iSM/bjsM and 
then Z^jY strongly depends on the superbubble radius. 
If so, the index, <5, in Equation [S] for helium is different 
from that for carbon, and hence the spectra of escaping 
CRs differ between helium and carbon. However, since 
the carbon to helium ratio in CRs provided by AMS-02 is 
constant within ^ 20 % in Figure [TJ the mixing fraction 
of ISM should be less than ^ 20 % at rigidity ii ^ 10 GV, 
and less than ^ 10 % at ii 10 TV. In other words, 
CRs originate primarily from the supernova ejecta and 
stellar winds in the interiors of superbubbles. 

For the Salpeter initial mass function [2^, massive 
stars with about 10 Mq mainly contribute to the CR pro¬ 
duction. According to a stellar nucleosynthesis model, 
stellar ejecta of a non-rotating massive star with 12 Mq 
give Zc^^-jYe^ = 0.0298 = 2.39 (^c,ism/Vism) [13, which 
is smaller than the value in the above argument. In this 
case, the mixing fraction of ISM may be about 50 %. 
However, we should keep in mind that the theoretical 
models of the stellar nucleosynthesis have still large un¬ 
certainties. 


V. ORIGIN OF SPECTRAL BREAKS 

Theoretical explanations for the spectral breaks in all 
the CRs at a single rigidity are classified into three types: 
the breaks produced intrinsically at the CR sources, the 
breaks produced during the propagation of CRs in the 
ISM, and the breaks due to the contribution of a local 
source at low/high rigidities. In the following, we will 
find that even the current AMS-02 data is not accurate 
enough to pin down the model for the breaks. 


A. Intrinsic 

There are several mechanisms to produce the breaks at 
the sources. As discussed in Section |H1 the Mach number 
dependence of the spectral index can produce the breaks 
within the framework of the standard DSA. The reac¬ 
celeration of pre-existing CRs at the superbubble [13 or 
the nonlinear acceleration of CRs [l3 lH, 113 Eire also 
the scenarios that give spectral breaks at the accelera¬ 
tion sites. Even when the CR spectrum at a source has 


a simple power-law form, if multiple kinds of sources [ 28 l | 
or sites [13 contribute to the CR flux or each spectral 
index has a dispersion [13, the superposed spectra can 
have breaks. B/C would not have a break in these ’’in¬ 
trinsic” scenarios, although there is an exception if the 
secondary boron is produced in the source (^ . 


B. Propagation 

If the energy dependence of the diffusion coefficient has 
a break at the rigidity ~ 300 GV, the observed CR spec¬ 
trum would have a break at the same rigidity [13 • In 
this case, the energy spectrum of B/C would also have 
a break, and a larger ratio is predicted in the higher en¬ 
ergy range Q- Although such a feature has not been 
confirmed by AMS-02, we cannot rule out this scenario 
because of the large uncertainty of data in the high en¬ 
ergy range. The situation is similar in other propagation 
models, e.g., a model in which the diffusion coefficient is 
not separable in energy and space [13 , and in which CRs 
are reaccelerated before escaping from our Galaxy [Hjj . 


C. Local source 

If there is a local CR source in addition to the Galactic 
sources, the observed CR spectrum would have a break 
at a certain rigidity, below/above which the local source 
dominates the CR flux [SJ]. The spectral difference be¬ 
tween the local source and the Galactic sources may be 
produced by such as the GR confinement [13 ■ However, 
0 show that the low-energy local source scenario fails 
to reproduce the observed B/G at low rigidities. On the 
other hand, in the high-energy local source scenario, B/C 
would be suppressed at high rigidities, where the error- 
bars of data by AMS-02 are still too large to rule out this 
scenario. In these scenarios, the ratio of CR protons to 
CR helium (and possibly to CR carbon) produced by the 
local source should coincide with that by the Galactic 
sources. 


VI. SUMMARY AND DISCUSSION 


We discussed that superbubbles are the most plausible 
source of Galactic CRs to explain the recent results of 
AMS-02, especially the helium and carbon hardenings. 
The spectra of CR helium and carbon are different from 
that of CR protons because the abundance profile is not 
uniform in superbubbles. All CR spectra may also break 
at i? ^ 300 GV since the shock Mach number evolves in 
high temperature 10® K) plasma. 

The superbubble origin of GRs is well motivated by 
other observations [22,(33. First, superbubbles contain 
^ 75% of all Galactic supernovae [^. Second, the ob¬ 
served abundances of light elements, Li, Be and B, in 
halo stars require CR carbon and oxygen in the early 
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metal-poor Galaxy and imply the superbubble origin [38l| . 
Third, the excess of the isotopic ratio ^^Ne/^°Ne sug¬ 
gests a link to winds from Wolf-Rayet stars, probably a 
major source of ^^Ne, and thus to superbubbles [s^ , al¬ 
though this is not supported by the updated wind model 
[13 • Finally, gamma rays are detected from OB asso¬ 
ciations such as in the Cygnus region [dlj and Wester- 
lund 2 [l^. Also theoretically, the particle acceleration 
is efficient in superbubbles . CR observations for 

^®Ni and ®°Fe suggest that CRs must be accelerated in 
10^ yr < t < 10^ yr after a supernova explosion (dl. l45l|. 
that is also consistent with the superbubble origin of CRs 

EF 


We suggested that the supernova ejecta dominate the 
ISM in the acceleration site to explain the AMS-02 obser¬ 
vations of He/p and C/p up to i? ^ 1 TV in Section lTVDI 
On the other hand, some studies based on CR data up 
to i? ~ 1 GV imply that the supernova ejecta fraction is 
~ 15-25% and the other is mixed from the ISM [l^.l4^. 
partly because this mixing ratio improves the power-law 
fit to the relation between the abundance and atomic 
mass. Because these arguments are based on the data in 
different rigidities, the mixing ratio derived in this paper 
{R > 30 GV) is allowed to be different from that derived 
in previous studies (i? < 1 GV). Thus, higher-energy 
CRs, observed by AMS-02, could be produced in more 
enriched core of superbubbles. 

If the composition of CRs with i? > 30 GV at the 
source is different from that with R < 1 GV, there must 
be some structures in CR spectra below i? = 30 GV. 
However, it is not always observable. As discussed in the 
last paragraph of the section lIVDi in order to explain 
the new AMS-02 data, the mixing fraction of supernova 
ejecta may be about 50 %. If so, the difference of the 


fraction between i? < 1 GV and i? > 30 GV is only by a 
factor of 2. Interestingly, data of AMS-02 actually sug¬ 
gest that the CR carbon to proton ratio, C/p, has a sig¬ 
nature of deviation from a single power law at i? ~ 10 GV 
(See Figure [T]). This is a kind of signature that supports 
the change of composition. In addition, structures in CR 
spectra below R < 30 GV could be smeared by the sec¬ 
ond order acceleration during the Galactic propagation. 
In Oder to discuss more quantitatively, we need to cal¬ 
culate the CR propagation in our Galaxy and take into 
account nonlinear effects, dust acceleration, solar modu¬ 
lation and so on [l^. These will be addressed in future 
work. 

In the near future, AMS-02 will provide precise data for 
heavy nuclei, which are predicted to have similar spectra 
in our model (possibly except for CR iron). Furthermore, 
the precise data for heavy nuclei with i? < 30 GV will 
give us some hints for the composition ch ang e between 
R < 1 GV and R > 30 GV. CALET [if and ISS- 
CREAM [i^ will also observe CRs above R ~ 1 TV. 
The B/C provided by future experiments can establish 
whether the breaks at R ~ 300 GV are due to intrinsic, 
propagation, or local source origin. 
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